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Abstract.  In  the  present  investigation,  slip-continuum  and  molecular  (DSMC)  flow  models  have  been  used  to  compare 
the  drag  coefficient  on  a  confined  microsphere  over  a  range  of  Knudsen  numbers  between  0.01  and  0.4.  The  simulations 
consider  a  fixed  Reynolds  number  (5e=0.125)  and  a  range  of  blockage  ratios  from  2.0  up  to  13.0.  Two  separate  problem 
formulations  have  been  studied;  the  first  considers  flow  past  a  confined  sphere  moving  with  a  constant  velocity  along  the 
central  axis  of  the  pipe  while  the  second  considers  Poiseuille  flow  past  a  stationary  confined  sphere.  Finally,  their 
combination  -  the  case  of  a  confined  sphere  moving  within  the  carrier  Poiseuille  flow  is  investigated  using  the  DSMC 
method. 


INTRODUCTION 

The  emergence  of  Micro-Electro-Mechanical-Systems  (MEMS)  as  a  key  enabling  technology  has  led  to  the 
development  of  an  increasing  number  of  gas-phase  microfluidic  systems.  Potential  applications  are  numerous  and 
include  miniaturised  heat  exchangers,  portable  gas  chromatography  systems,  miniaturised  gas  sensors  and  novel 
high-throughput  gas  flow  cytometers.  Elowever,  one  of  the  most  important  issues  influencing  MEMS  research  is  the 
growing  realisation  that  gas-phase  microflows  are  dominated  by  non-continuum  or  rarefaction  effects. 

Previous  numerical  studies  by  the  authors  have  investigated  the  important  problem  of  low  Reynolds  number  slip 
flow  past  a  confined  microsphere  within  a  circular  pipe  [1,2].  The  simulations  were  conducted  using  two  different 
numerical  approaches  based  upon  continuum  (Navier-Stokes)  [1]  and  molecular  (DSMC)  [2]  model  descriptions. 
The  studies  have  resulted  in  a  number  of  fundamental  questions  regarding  the  application  of  both  conventional 
(continuum)  and  molecular  (DSMC)  models  for  gas  microflows.  One  such  question  concerns  the  upper  limit  of 
Knudsen  number  for  which  a  continuum  slip-boundary  solution  might  be  valid.  Another  is  whether  the  DSMC 
method  is  able  to  cover  a  sufficiently  wide  range  of  Knudsen  numbers  so  that  continuum  and  molecular  data  can  be 
compared. 

In  the  present  investigation,  slip-continuum  and  molecular  models  have  been  used  to  compare  the  drag 
coefficient  on  a  confined  microsphere  over  a  range  of  Knudsen  numbers  between  0.01  and  0.4.  The  simulations 
consider  a  fixed  Reynolds  number  (Re=0.  1 25 )  and  a  range  of  blockage  ratios,  B  =  (pipe  diameter)/( microsphere 
diameter)  from  5=2.0  up  to  5=13.0.  Two  separate  problem  formulations  have  been  studied;  the  first  considers  flow 
past  a  confined  sphere  moving  with  a  constant  velocity  along  the  central  axis  of  the  pipe  while  the  second  considers 
Poiseuille  flow  past  a  stationary  confined  sphere.  The  continuum  model  employs  a  tangential  slip-velocity  boundary 
condition  and  assumes  that  compressibility  effects  can  be  neglected  while  the  molecular  simulations  are  conducted 
using  a  standard  DSMC  method.  The  DSMC  approach  allows  non-equilibrium  and  rarefaction  effects  to  be  captured 
over  a  wide  range  of  Knudsen  numbers  but  has  the  disadvantage  of  requiring  extremely  large  numbers  of  particles  at 
low  Knudsen  numbers,  making  the  technique  prohibitively  expensive  in  terms  of  computational  cost.  The  reason  is 
that  the  macroscopic  gradients  and  bulk  velocity  in  micro-flows  of  practical  interest  are  usually  small  and  the 
statistical  error  of  the  DSMC  calculations  is  comparable  to  the  magnitude  of  the  macroscopic  flow  characteristics. 
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The  DSMC  calculations  have  shown  that  these  difficulties  can  be  overcome,  under  certain  conditions,  by  using  a 
simple  data  filter  for  computing  the  averaged  flow  fields.  It  should  be  noted  that  a  promising  alternative  approach  is 
the  information  preservation  method  proposed  by  Fan  and  Shen  [3],  In  the  present  paper,  the  DSMC  approach  has 
been  used  to  simulate  flows  down  to  a  Knudsen  number,  Kn  =  (mean  free  path)/( sphere  diameter)  of  0.01  for  a 
blockage  ratio  of  5=2.0,  Kn= 0.02  for  5=5.0,  and  Kn=0.04  for  5=13.0  (approaching  the  case  of  an  unconfined 
sphere).  For  this  last  case,  the  drag  force  has  also  been  compared  against  the  analytical  expression  for  the  drag  on  an 
unconfined  sphere  presented  by  Beresnev  et  al.  [4].  In  this  case,  the  continuum  slip-flow  calculations  overpredict  the 
drag  coefficient  at  all  Knudsen  numbers  but  converge  to  the  Beresnev  solution  as  Kn  tends  to  zero.  The  DSMC 
results  show  good  agreement  with  Beresnev  et  a/.'s  analytical  solution  for  Kn  >  0.04.  In  the  case  of  a  confined 
sphere,  the  continuum  model  predicts  a  larger  drag  coefficient  than  the  DSMC  method  when  Kn  is  larger  than  some 
critical  value  that  depends  on  the  blockage  ratio.  However,  at  very  small  Knudsen  numbers,  the  DSMC  simulations 
predict  a  larger  drag  coefficient  than  the  continuum  slip  model.  This  interesting  result  may  be  explained  by  the 
effects  of  gas  compression  in  front  of  the  microsphere  and  rarefaction  behind  it. 

Finally,  we  report  some  preliminary  DSMC  results  concerning  a  more  realistic  case  of  flow  past  a  confined 
sphere  moving  within  the  carrier  fluid  along  the  central  axis  of  the  pipe.  In  other  words,  we  have  generalized  the 
problem  and  considered  the  combination  of  the  two  limiting  formulations  mentioned  above.  The  combined 
formulation  makes  it  possible  to  estimate  the  relative  motion  between  the  sphere  and  the  bulk  flow. 


PROBLEM  FORMULATION 

We  consider  a  monatomic  gas  flow  past  a  confined  microsphere  with  diameter,  D,  in  a  circular  pipe  with 
diameter,  H  (H>  D).  The  flow  is  simulated  using  an  axially-symmetric  computational  domain  with  a  finite  length,  L. 
The  computational  domain  and  the  imposed  local  cylindrical  coordinate  system  (x,r)  are  shown  in  Fig.  1.  Initially, 
the  pipe  is  filled  with  gas  of  density,  p0=mn0  (m  is  the  mass  of  the  gas  molecules  and  n0  is  the  number  density)  and 
temperature,  T0.  The  pipe  wall  and  the  inlet/outlet  gas  fluxes  are  assumed  to  have  the  same  temperature,  T0.  Diffuse 
reflection,  or  a  corresponding  slip-velocity  boundary  condition,  is  assumed  on  both  the  pipe  and  particle  surfaces. 
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FIGURE  1.  Local  computational  domain  around  the  sphere.  Inflow  velocity  profile  Uo  and  boundary  velocity  C/w  for:  case  (a),  a 
confined  sphere  moving  with  a  constant  velocity  along  the  central  axis;  case  (b),  Poiseuille  flow  past  a  stationary  confined  sphere; 
case  (c),  a  confined  sphere  moving  with  velocity  U(t)  in  the  Poiseuille  flow. 


Under  these  general  conditions  we  have  studied  three  flow  regimes.  The  first  case  concerns  a  confined  sphere 
moving  with  a  constant  velocity  U0  along  the  central  axis  (as  shown  in  Fig.  1(a)).  In  the  local  coordinate  system  this 
means  that  the  sphere  is  at  rest  while  the  pipe  wall  is  moving  with  velocity  C/w=  U<,.  The  inflow  at  x=0  has  a  uniform 
velocity  profile,  U0.  At  the  outlet,  the  axial  gradients  of  the  velocity  components  are  assumed  to  be  zero. 

The  second  case  concerns  Poiseuille  flow  past  a  stationary  confined  sphere  (Fig.  1(b)).  In  the  local  coordinate 
system,  the  wall  velocity  Uw=0.  At  the  inflow  boundary  (x=0),  a  fully-developed  slip-velocity  profile  [1]  is 
prescribed  parallel  to  the  longitudinal  axis  of  the  pipe: 
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where  U  is  the  mean  velocity  in  the  pipe  and  a  is  the  tangential  momentum  accommodation  coefficient  (a  =  1  for 
diffuse  reflection).  At  the  outlet,  the  axial  gradients  of  the  velocity  components  were  assumed  to  be  zero. 

The  third  case  (Fig.  1(c))  is  a  combination  of  the  first  two  and  describes  the  microsphere  motion  along  the  central 
axis.  The  motion  is  induced  by  the  drag  of  the  carrier  Poiseuille  flow.  The  instantaneous  velocity  of  the  sphere  U(t) 
is  computed  from  the  equation: 

mpdL4f^  =  FrAt\  (2) 

at 

where  mp  is  mass  of  the  sphere  and  FD(t)  is  the  instantaneous  drag  force  acting  on  the  sphere.  For  simplicity,  in  Eq. 
(2)  we  have  neglected  Brownian  motion  and  assumed  there  are  no  external  forces  on  the  sphere.  The  initial  condition 
is  {/( 0)=0.  In  the  imposed  coordinate  system,  the  inflow  condition,  Eq.  (1),  at  x=0  is  modified  to  give 

U0(t)  =  U0  —U(t).  (3) 


The  wall  velocity  thus  becomes  Uw  =  -U(t).  For  clarity  in  the  computational  algorithm  used  for  case  (c),  it  should  be 
noted  that  the  longitudinal  component  of  the  velocity  field,  U{trX,r),  computed  in  the  local  coordinate  system,  was 
also  modified  each  time  step: 

U(t  +  At,  x,  r)  =  U(t,  x,  r )  —  ^t.  (4) 

dt 


Numerical  simulations  were  used  to  assess  the  drag  experienced  by  a  confined  sphere  exposed  to  low  Reynolds 
number  flow.  In  case  (a),  the  Reynolds  number.  Re,  was  defined  using  the  constant  velocity,  U0,  and  the  sphere 
radius,  a=D/2,  as  the  velocity  and  length  scales: 


Re  = 


paU  o 

9 

A 


(5) 


where  fi  is  the  viscosity  of  the  gas,  whilst  in  cases  (b)  and  (c),  Re  was  defined  using  the  mean  velocity  in  the  pipe, 
U  ,  instead  of  U0.  The  second  parameter,  the  Knudsen  number,  Kn=X/D,  determines  the  degree  of  rarefaction  of  the 
gas.  The  mean  free  path  of  molecules,  X,  was  defined  for  an  ideal  gas  modeled  as  rigid  spheres  with  diameter  d  (Eq. 
(9)).  Blockage  effects  were  studied  by  varying  the  ratio  between  the  diameter  of  the  confining  pipe,  H,  and  the 
diameter  of  the  sphere,  D. 


COMPUTATIONAL  CONSIDERATIONS 

For  flow  cases  (a)  and  (b),  slip-continuum  and  molecular  models  have  been  used  to  compare  the  drag  coefficient 
on  a  confined  microsphere  over  a  range  of  Knudsen  numbers  between  0.01  and  0.4.  All  simulations  consider  a  fixed 
Reynolds  number  (Re=0.  1 25)  and  a  range  of  blockage  ratios  from  B= 2.0  up  to  .8=13.0  (for  the  molecular  model)  and 
from  .8=2.0  up  to  8=40.0  (for  the  slip-continuum  approach).  The  continuum  model  assumes  that  compressibility 
effects  can  be  neglected  and  employs  a  tangential  slip-velocity  boundary  condition  at  the  pipe  wall  and  surface  of 
the  microsphere: 

2-a  KnD  ((,) 

u,  = - r,  ’ 

o  p 

where  rt  is  the  shear  stress  at  the  surface.  The  governing  hydrodynamic  equations  [1]  were  solved  using  a  two- 
dimensional  finite-volume  Navier-Stokes  solver  developed  by  the  Computational  Engineering  Group  at  CCLRC 
Daresbury  Laboratory. 

The  molecular  simulations  were  conducted  in  a  two-dimensional  cylindrical  domain  (x,r)  e  (0,/.)x(0,///2)  using 
standard  DSMC  procedures  for  simulation  of  flows  with  axial  symmetry  [5].  The  computational  domain  was 
covered  by  a  uniform  rectangular  basic  grid  with  NxxNr  cells.  For  the  most  difficult  computational  case  (Kn= 0.04, 
8—  1 3.0),  a  pipe  domain  with  dimensions  (T,R)=( 64x20)  pm  was  employed  and  the  radius  of  the  spherical  particle 
was  taken  as  Rp=1.55  pm.  The  total  number  of  simulated  molecules  in  the  pipe  volume  was  approximately  20xl06. 
In  this  case,  the  computational  domain  was  covered  by  a  grid  with  2400x800  basic  cells.  The  basic  cells  near  the 
surface  of  the  sphere  were  subdivided  dynamically  into  subcells  in  order  to  meet  the  spatial  resolution  requirements 
of  the  method.  Flow  field  sampling  was  conducted  on  a  coarser  grid  with  larger  cells  containing  20x20  basic  cells. 
This  multilevel  grid  scheme  allowed  the  calculation  of  the  molecular  processes  correctly  on  an  adaptive  fine  grid  and 
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at  the  same  time  provided  a  meaningful  sample  size  for  the  macroscopic  variables.  In  addition,  a  simple  filtering 
procedure  was  applied  along  with  standard  time  averaging  of  the  accumulated  flow  field  data.  For  all  cases,  the 
computed  temperature  was  7,0=288  K.  The  smallest  time  step  used  for  computing  the  near-continuum  cases  was 
At=1.0xl0  10  s.  The  force  acting  on  the  sphere  was  determined  from  the  rates  of  delivery  and  removal  of  momentum 
by  incident  and  reflected  molecules  from  the  particle  surface  within  a  time  interval,  A ts: 

<  F  >=  —  T  m(r<+>  -  Vj-> )  •  (V 

A  ts  ^  V  *  ’ 

In  case  (c),  the  instantaneous  value  of  the  drag  force  was  taken  to  be  FD(i)  =<Fx(t)>  estimated  over  a  small 
macroscopic  scale  time  interval  A4=(?-ts)=(l~5)102A8  In  cases  (a)  and  (b),  the  drag  force  was  averaged  over  a  large 
time  interval,  Ats=  I  O'Al,  after  reaching  the  steady-state  flow  conditions. 


NUMERICAL  RESULTS 


The  slip-continuum  Navier-Stokes  (NS)  and  molecular  model  (DSMC)  results  were  compared  for  two  separate 
flow  problems  (case  (a)  and  case  (b)).  The  Knudsen  number,  Kn,  and  the  blockage  ratio,  B,  were  varied  while  the 
Reynolds  number  was  fixed  at  a  small  value,  Re=0.125.  Since  the  slip-continuum  model  assumes  a  constant  density, 
the  variation  of  the  Knudsen  number  does  not  affect  the  Reynolds  number  and,  consequently,  for  fixed  Re,  there  is 
no  correlation  between  gas  density  and  flow  velocity  when  changing  Kn.  This  is  not  the  case  with  the  molecular 
model.  All  DSMC  results  were  obtained  in  dimensional  form  using  the  “hard  sphere”  model  with  a  molecular  mass 
of  4.815xl0'26  kg  and  a  molecular  diameter  of  3.7xlO"10  m.  Under  these  conditions,  the  mean  free  path  variation  for 
a  gas  in  equilibrium  depends  only  on  the  gas  density.  Thus,  to  keep  Re  constant,  when  changing  Kn,  the  product 
(pU)  must  be  also  kept  constant.  It  can  easily  be  shown  that  the  following  relation  between  Re,  Kn  and  Mach 
number,  Ma,  is  valid  for  a  monatomic  hard  sphere  gas: 


Re  = 


I  32  Ma 
V  1 5zr  Kn 


(8) 


where  the  Mach  number  and  the  Chapman-Enskog  first  approximation  for  viscosity  and  mean  free  path  are 


Ma  = 
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V(5/3  )kTJ  m 


/u  =  — A 
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2rtkT., 


A  =  (V 2^ii2;i0 )  1  • 
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Equation  (8)  indicates  that  to  keep  Re  constant  while  changing  Ain,  necessitates  changing  the  Mach  number,  Ma. 


Case  (a):  Confined  Sphere  Moving  with  a  Constant  Velocity 

For  case  (a),  the  slip-continuum  Navier-Stokes  (NS)  calculations  were  performed  for  a  range  of  Knudsen 
numbers  from  Kn= 0.0  (continuum  flow  regime)  to  Kn= 0.4  and  a  range  of  blockage  ratios  from  FFD= 2  to  H/D= 40. 
The  latter  can  be  considered  as  approaching  the  case  of  an  unconfined  sphere.  The  molecular  (DSMC)  approach  was 
used  to  simulate  flows  down  to  a  Knudsen  number,  Kn= 0.01  for  a  blockage  ratio  of  B= 2.0,  Kn= 0.02  for  .8=5.0,  and 
A7z=0.04  for  IF  1 3.0  (approaching  the  case  of  an  unconfined  sphere).  We  compare  the  normalized  total  drag 
coefficient  on  the  microsphere,  CoRe/12  [1]  versus  the  Knudsen  number  for  three  separate  blockage  ratios,  8=2.0, 
.8=5.0  and  IF  1 3.0  (DSMC)  and  8—40.0  (NS).  For  the  largest  value  of  8,  the  normalized  drag  coefficient  was  also 
compared  against  the  analytical  formula  presented  by  Beresnev  et  al.  [4],  which  is  in  good  agreement  with  other 
theoretical  and  numerical  results  (see  for  example  [6,7]).  For  complete  diffusive  reflection  on  the  surface  of  the 
sphere,  the  normalized  analytical  drag  coefficient  [4]  can  be  written  as  follows: 

CD  Re  _  8  +  n 
12  ~~  18 

The  results  of  the  comparison  are  shown  in  Fig.  2(a).  The  DSMC  calculations  were  extended  up  to  AT«=10.0.  For 
B»  1  (approaching  the  case  of  an  unconfined  sphere)  the  DSMC  results  were  in  good  agreement  with  Beresnev’s 
formula  (10).  An  exception  is  the  data  point  at  Kn=0.QA  where  a  meaningful  sample  size  of  DSMC  data  could  not  be 
reached.  The  reason  for  the  observed  difference  at  large  Kn  is  the  formation  of  a  shock  wave  when  the  bulk  velocity 
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exceeds  the  sonic  speed.  In  the  case  of  an  unconfined  sphere,  the  continuum  slip-flow  calculations  overpredict  the 
drag  coefficient  at  all  Knudsen  numbers  but  converge  to  the  Beresnev  solution  as  Kn  tends  to  zero.  In  the  case  of  a 
confined  sphere  (5=5.0  and  5=2.0),  the  continuum  model  predicts  a  larger  drag  coefficient  than  the  DSMC  method 
when  Kn  is  larger  than  some  critical  value  that  depends  on  the  blockage  ratio.  However,  at  very  small  Knudsen 
numbers,  the  DSMC  simulations  predict  a  larger  drag  coefficient  than  the  continuum-slip  model.  This  interesting 
result  may  be  explained  by  the  effects  (observed  in  the  DSMC  data.  Fig.  2(b))  of  gas  compression  in  front  of  the 
microsphere  and  rarefaction  behind  it. 


FIGURE  2.  (a)  Normalized  drag  coefficient  versus  Knudsen  number:  slip-continuum  model  (dashed  lines);  DSMC  data  (circles). 
For  5»1  both  slip-continuum  and  DSMC  results  are  compared  with  the  analytical  formula,  Eq.  (10)  (solid  line);  (b)  Density 
distribution  along  the  x-axis:  DSMC  data  (circles),  averaged  density  in  front  of  the  sphere  (solid  line),  averaged  density  behind 
the  sphere  (dashed  line). 


Case  (b):  Gas  Flow  Past  a  Stationary  Confined  Sphere 

The  normalized  drag  coefficient  obtained  for  Poiseuille  flow  past  a  stationary  confined  sphere  (Fig.  1(b))  is 
shown  in  Fig.  3(a).  The  Knudsen  number  in  the  DSMC  simulations  ranged  from  0.01  to  0.6  while  the  slip- 
continuum  model  calculations  were  performed  over  the  ranges  shown  in  Fig.  3(a).  Two  aspect  ratios,  5=2.0  and 
5=5.0,  were  considered.  The  results  for  5=2.0,  case  (b)  (solid  lines),  are  given  together  with  the  corresponding 
results  for  case  (a)  (dash-dot  lines).  Again,  as  in  case  (a),  the  DSMC  computed  normalized  drag  coefficient  is  larger 
than  the  corresponding  slip-continuum  model  at  small  Knudsen  numbers,  and  less  than  the  continuum  approach  at 
large  Knudsen  numbers.  An  interesting  result  observed  in  the  DSMC  computed  drag  coefficient  is  the  larger  slope  of 
the  curve  in  case  (b)  compared  to  case  (a).  In  contrast,  the  slope  of  the  slip-continuum  results  is  almost  the  same  in 
both  cases.  This  difference  is  demonstrated  in  Fig.  3(b)  which  shows  the  dependence  of  the  dimensional  drag  force 
on  the  Knudsen  number  for  Re=0.  1 25  and  5=2.0.  It  can  be  seen  that  the  different  flow  gradients  around  the  sphere 
have  the  opposite  effect  on  the  drag  force  variation.  In  case  (a),  the  drag  is  approximately  proportional  to  Kn  while 
in  case  (b),  the  drag  is  inversely  proportional  to  Kn. 

Case  (c):  Confined  Sphere  Moving  with  Velocity  U(t)  in  the  Carrier  Poiseuille  Flow 

Figure  4(a)  illustrates  the  time  evolution  of  the  particle  velocity  in  the  pipe  under  the  action  of  the  carrier 
Poiseuille  flow  for  case  (c).  The  results  were  obtained  for  Kn=0A  using  the  DSMC  method.  From  Fig.  4(a)  one  can 
conclude  that  the  blockage  ratio,  5,  strongly  affects  the  particle  velocity.  The  particle  velocity  was  computed  at  each 
time  step  from  Eq.  (2),  where  the  mass,  mp,  was  taken  for  a  sphere  with  density  1800  kg/m3.  The  drag  force  (Eq.  (7)) 
oscillations  (circles)  are  illustrated  in  Fig.  4(b)  for  5=5.0.  When  the  filtered  mean  drag  (thick  solid  line)  tends  to 
zero,  the  corresponding  velocity  asymptotically  reaches  its  constant  value. 
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(a)  (b) 

FIGURE  3.  (a)  Comparison  of  normalized  drag  coefficient  versus  Knudsen  number  for  the  DSMC  data  (triangles)  and  the  slip- 
continuum  model  (lines)  for  case  (a)  (dash-dot  lines)  and  case  (b)  (solid  lines);  (b)  Dimensional  drag  versus  Kn,  for  the  DSMC 
data  (triangles)  for  both  case  (a)  and  case  (b). 


(a) 


(b) 


FIGURE  4.  (a)  Time  evolution  of  velocity  of  a  confined  microsphere  dragged  by  the  gas  flow.  The  results  are  shown  for  two 
blockage  ratios  B= 2.0  (circles)  and  3=5.0  (triangles);  (b)  Drag  oscillations  (circles)  and  filtered  mean  drag  (solid  line)  for  5=5.0. 


CONCLUSIONS 

Slip-continuum  and  molecular  models  have  been  used  to  compare  the  drag  coefficient  on  a  confined  microsphere 
over  a  range  of  Knudsen  numbers  between  0.01  and  0.4.  The  results  show  that  the  continuum  model  predicts  a  larger 
drag  coefficient  than  the  DSMC  method  when  Kn  is  larger  than  some  critical  value  that  depends  on  the  blockage 
ratio.  However,  for  small  Kn,  the  DSMC  predictions  indicate  a  larger  drag  coefficient  than  the  continuum  model. 
This  interesting  result  may  be  explained  by  the  effects  of  gas  compression  in  front  of  the  microsphere. 


REFERENCES 

1.  Barber,  R.  W.,  and  Emerson,  D.  R.,  In  Rarefied  Gas  Dynamics:  23rd  Int.  Symp.-2002,  edited  by  A.D.  Ketsdever  and  E.P. 
Muntz,  A1P  Conference  Proceedings  663,  American  Institute  of  Physics,  New  York,  2003,  pp.  808-815. 

2.  Ota,  M.,  Kuwahara,  K.,  and  Stefanov,  S.,  In  Rarefied  Gas  Dynamics:  23rd  Int.  Symp.-2002,  edited  by  A.D.  Ketsdever  and 
E.P.  Muntz,  AIP  Conference  Proceedings  663,  American  Institute  of  Physics,  New  York,  2003,  pp.  776-783. 

3.  Fan,  J.,  and  Shen,  C.,  J.  Comp.  Phys.  167,  393-412  (2001). 

4.  Beresnev,  S.  A.,  Chernyak,  V.  G.,  and  Fomyagin,  G.  A.,  J.  Fluid  Mechanics  219,  405-421  (1990). 

5.  Bird,  G.  A.,  Molecular  Gas  Dynamics  and  the  Direct  Simulation  of  Gas  Flows,  Clarendon  Press,  Oxford,  1994. 

6.  Cercignani,  C.,  Pagani,  C.D.,  and  Bassanini,  P.,  Phys.  Fluids  11,  1399-1403  (1968). 

7.  Sone,  Y.,  and  Aoki,  K.,  “Forces  on  a  Spherical  Particle  in  a  Slightly  Rarefied  Gas”,  In  Rarefied  Gas  Dynamics:  1 0th  Int. 
Symp.  —  1976,  edited  by  J.  Potter,  AIAA  Progress  in  Astronautics  and  Aeronautics  51,  New  York,  1977,  pp.  417-433. 


706 


